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Abstract 

An Ising model for damage spreading with a probability of damage healing {q = 1 ~ p) is proposed and 
studied by means of Monte Carlo simulations. In the limit p 1 the new model is mapped to the standard Ising 
model. It is found that, for temperatures above the Onsager critical temperature (Tc), there exist a no trivial 
finite value of p that sets the critical point (pc) for the onset of damage spreading. It is found that Pc depends 
on T, defining a critical curve at the border between damage spreading and damage healing. Transitions along 
such curve are found to belong to the universality class of directed percolation. The phase diagram of the model 
is also evaluated showing that for large T one has Pc oc {T — Tc)°', with a = 1. Within the phase where the 
damage remains active, the stationary value of the damage depends lineally on both p — Pc and T — Tc- 
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I. Introduction 



The damage spreading method has been widely used to study the critical properties of Ising-like systems (TJ [21 [31 
HI HI [SI [ll [8] as well as to spin glass [9l[T0]. The method is based on the synchronous Monte Carlo update of two 
distinct spin configurations that are evolving from an almost identical initial state [Tl] . Since the Ising model lacks 
of an intrinsic dynamic, it has to be chosen a particular one. Among others, the more frequently used are Glauber 
and heat-bath dynamics. In contrast to usual statistical Monte Carlo studies, damage spreading results depends 
on the used dynamic pj] . This behavior has been used to evaluate the dynamical exponent z of the 2D and 3D 
Ising models with different dynamics |3l[T2l[l3]. Also, very recently, it has been shown that damage spreading is 
a powerful and useful technique for the numerical study of the role of the interfaces between magnetic domains on 
the propagation of perturbations in magnetic materials [HI [15]. Within this context, we have reported that the 
presence of interfaces act as a "catalyst" of the damage in at least two different ways: speeding up the propagation 
and causing an enhancement of the total damaged area [HI [15] . 

The basic idea for the implementation of the damage spreading method [11] is to start from an equilibrium 
configuration of the system at temperature T, which is generically called iS'^(r). Subsequently, a very small 
perturbation is applied to that configuration in order to obtain a new one, i.e. the so called perturbed configuration 
S^{T). Usually the perturbed configuration is obtained by flipping a small number of spins of the unperturbed 
conflguration. Then one has to study the time evolution of the perturbation in order to investigate under which 
conditions such a small perturbation will grow up indefinitely or eventually (hopefully) it will vanish and become 
healed. In order to follow the time evolution of the perturbation an useful method is to measure the "Hamming 
distance" or damage between the unperturbed and the perturbed configurations. The total damage D(t) is defined 
as the fraction of spins with different orientations, that is 



where the summation runs over the total number of spins N and the index I {1 < I < N) is the label that identifies 
the spins of the configurations. Starting from a vanishing small perturbation D{t = 0) ^ one can expects at 
least two scenarios, namely: a) D{t oo) and the perturbation is irrelevant because it become healed and 
b) D{t oo) goes to some well defined finite value. In the latter, frequently undesired case, the perturbation is 
relevant because it can not be healed out. For further details see the review |11) . 

The aim of this work is propose and study a model for damage spreading with damage healing, which is based in 
the two dimensional {2D) Ising model. It is well known that the 2D Ising magnet undergoes a second-order order- 
disorder transition when the temperature is raised from a relatively low initial value. The location of the critical 
point is know exactly and it is given by the so called Onsager critical temperature — 2.269.... According to 
previous studies of damage spreading using Glauber dynamics it is known that there is a critical damage temperature 
given by Td = 0.992 Tq [T6lll7j. such as for T > Tjj the damage spreads out over the whole sample while for T < Tjj 
the damage becomes healed after some finite time. Consequently, the proposed model, that incorporates a healing 
probability, is suitable for the study of damage spreading above Td , as well as to gain insight on the robustness of 
the damage behavior in the Ising model. 

The manuscript is organized as follows: in Section II we propose the model and describe the numerical procedure 
for the simulation of damage spreading. Section III is devoted to the presentation and discussion of the results, 
while our conclusions are stated in Section IV. 

II. A Model for Damage Spreading with Healing and the Monte Carlo 
Simulation Method. 

The model is based on the 2D ferromagnetic Ising model, which for a square lattice of side L can be described by 
the following Hamiltonian H: 




(1) 




(2) 
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where cr^ is the spin variable, corresponding to the site of coordinates (i, j), that may assume two different values, 
namely = ±1, J > is the coupling constant of the ferromagnet and the summation of ((2]) runs over all the 
nearest-neighbor pairs of spins. 

For the purpose of the simulations, periodic boundary conditions are always used. Furthermore, we have used 
the Glauber dynamics. In order to implement this dynamic a randomly selected spin is flipped with probability 
Piflip) given by: 

exp(~/? ■ AH) 
^^^'^^^ = l + cxp(-/3.Ag) 

where AH is the difference between the energy of the would-be new configuration and the old configuration, and 
f3 = 1/kT is the usual Boltzmann factor. The temperature is measured in units of the Onsager critical temperature 
of the 2D Ising model. 

In order to study the time evolution of the damage spreading a meaningful definition of the Monte Carlo time 
step (mcs) is necessary. For this purpose the standard definition is adopted according to that during one mcs all 
L X L spins of the sample are flipped once, in the average. 

For the practical implementation of equation ([1]), first an equilibrium configuration (say configuration Sf) is 
generated. For this purpose one starts from an initial configuration with all TV = L x L spins oriented at random. 
The application of the Glauber dynamics leads to the desired equilibrium configuration after lO^mcs. Subsequently, 
a fully damaged replica of such configuration is created, such as is the mirror image of and consequently 
D{0) = 1. Of course, due to the spin fiip symmetry, the replica configuration is also equilibrated. 

In standard damage spreading studies, the standard Monte Carlo procedure is then applied to both configura- 
tions, where the same sites are randomly selected and the same random numbers are used in both systems in order 
to perform the updates. For the purpose of the proposed model we first define the healing probability q and we 
used the same method but according the following rules: 

i) A sample site of coordinates is selected at random. Then, 

iia) If S^{i,j) — S^{i,j), which means that the selected site is not damaged, the Glauber dynamics is followed 
according to the standard procedure [H]. However, 

iib) If S^{i,j) ^ S^{i,j), i.e. for a damaged site, a new random number h is generated and one proceeds as 
follows: 

iibl) If h < q, the damage is healed, so that one sets S^{i,j) — S^{i,j). However, 
iib2. li h> q, the standard dynamics is applied. 

According to the defined rules and taking p = 1 — q, one has that for p — > 0, the damage would become healed 
for all temperatures, while in the limit p — > 1, the usual damage spreading dynamics in the Ising model is recovered. 

Using the above described procedure, we have followed the time evolution of the damage D{t), which is evaluated 
according to equation (U), for different values of p, 1.1 < T < 50 and the lattice size 50 < i < 1000. Simulations 
are stopped when the damage is healed, otherwise they are performed up to i = 65000 mcs. 



III. Results and discussion 

Figure 1 shows log- log plots of D{t) versus t obtained at T = 2 for different values of p. It is found that: i) for p > Pc 
the damage becomes healed for finite times and the log-log plots of D{t) versus t exhibit a downward curvature, ii) 
For p < Pc the damage quickly propagates and the log-log plots of D{t) versus t show that the damage reaches a 
stationary value. Finally, just at p = Pc = 0.1895(5) the damage decays according to a power law. So, the change 
of the curvature of the plots shown in figure 1 allows us to identify the critical probability Pc for damage spreading. 
The straight line observed in figure 1 for pc suggests a power-law behavior given by 

D{t) oc f-^ (4) 

where S is an exponent. The best fit of the data gives S = 0.450 ± 0.005. Similar plots obtained for T > Tc (not 
shown here for the sake of space) exhibit the same behavior as that observed in figure 1 and the typical slope, 
averaged within the range 2 < T < 50, is given hy S = 0.45 ± 0.01. 

On the other hand, for p > Pc figure 1 shows a fast decay of the damage. As suggested by the study of spreading 
behavior [18], we proposed an exponential decay given by 
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D{t) oc exp{-t/T) 



(5) 



where r is the characteristic time of decay. 

In order to check the scaHng Ansatz given by equation |[5|), figure 2 shows log-Hneal plots of D{t) versus t 
obtained for p > Pc- The obtained straight lines for the long time behavior allow us to evaluate the dependence of 
r on p. 

It is found that r increases when approaching the critical probability. Furthermore, plots performed at different 
temperatures (no shown here for the sake of space) also exhibit the same behavior, independent of T. Again, using 
the spreading formalism [18], the following power-law can be proposed 

T (X (p - Pc)-"^^ (6) 

where is the correlation length exponent for damage propagation along the time direction. Figure 3 shows log-log 
plots of T versus p — Pc obtained at different temperatures. The typical averaged value resulting from the fits is 
given by v\\ = 1.25(5). 

It is worth mentioning that our estimations for the exponents, given hy 5 = 0.45(1) and = 1.25(5) are in 
agreement with the well known critical exponents of the universality class of Directed Percolation (DP), namely 
5DP = 0.4505(10) [Hj and v^^ = 1.295(6) [20j. So, these results strongly suggest that the the proposed model of 
damage spreading with damage healing belongs to the DP universality class. This statement is further supported 
by the fact that the limit p ^ 1 corresponds to the well known Ising model where the damage spreading transition 
is known to lie within the DP universality class [TBI [17]. 

On the other hand, as follows from figure 1, for p < p^ and after a short decay period {t « 10^ mcs), the damage 
becomes stabilized (or saturated) with typical average values Dsat that depend on p. This trend has been checked 
by means of extensive simulations run up to i = 10^ mcs. Furthermore, we have also checked that the stationary 
values of Dgat exhibit negligible finite-size effects for L > 200. Figure [4] shows the dependence of Dgat on {p^ — p) 
in a lineal-lineal plot. It follows that D^at increases lineally as a function of {p^ — p), and that the slope F{T) 
depends on temperature. This assumption is valid only in the limit oi p ^ p^ . Accordingly, we propose 

Dsat^Do + FiT){p^ -p), (7) 

where Dq « 0.05 is almost a constant independent of T (see figure |4]). The slope of the lines can be obtained by 
fitting the data of figure [4] and it is found that F(T) also depends lineally with T (see figure [5]) according to 

F{T) = Fo + a(^^^^ (8) 

where Fq w and A w 1.5 are also constants. 

So, by inserting equation ^ into equation ([7]), one obtains 



D,at = Do+(^F„ + A^^^y{p-pf), (9) 

that gives the full dependence of Dsat on both p and T. The scaling plot suggested by equation ^ holds acceptably 
as shown in figure El 

Performing plots of the time evolution of the damage using different values of p (as shown in figure [T]) and 
changing T, it is finally possible to evaluate the phase diagram of the model, i.e. a plot of p^ as a function of 
(T — Tc)/Tc, as shown in figure[7l It is found that p^ for very large values of T, as expected from the definition 
of the model. For the asymptotic decay of we propose 

pf-f^V" (10) 



where the best fit of the data, within the range T ^ 'V-'iiTc ~ 3.4 J, gives a = 1 for the exponent of equation l|10p . It 
is worth mentioning that, for the lattices used in the simulations, the phase diagram exhibits negligible deviations 
due to finite size effects, as shown in figure [7l On the other hand, close to Tc, it is observed that p^ reaches a 
maximum value of the order of p^ ~ 0.22 for T = 1.5rc ~ 3.4J. We have carefully checked that this behavior. 
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which remains to be understood, is not due to the operation of finite size effects. Finally, let us remark that all 
phase transitions at the critical curve at the border between damage healing and damage spreading, shown in figure 
[ll belong to the universality class of directed percolation. 

IV. Conclusions 

Based on the fact that in the 2D Ising model the damage is not healed above the critical temperature, we proposed 
a model that introduces a healing probability. It is found that for T > Tc the damage becomes healed for critical 
values of p. The phase diagram of the model is obtained and it is shown that the transitions between states of 
damage spreading and heahng belong to the universality class of directed percolation in (2 + 1)— dimensions. Our 
results thus support the conjecture that damage spreading transitions are generically in the universality class of 
directed percolation |T6], provided the fact that such transitions do not coincide with the critical point of the 
undamaged system. 
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Figure 1: Log-log plots of the damage {D{t)) versus time t. Data obtained using lattices of side L = 400, for 
T — 2.0Tc and taking different values of p, as listed in the figure. The full line has slope —S = —0.45 (see equation 
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Figure 2: Log- lineal plots of the damage {D{t)) versus t for data corresponding to p > pj? and obtained taken 
L = 400 and T = 2.0Tc. 



7 



10 



10 



10 



1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 


' ' 1 


1 1 1 


- 




• 


T=1.5 - 


- 


5 \ 

\ 


■ 


T=2.0 


- 


\ 

\I 


♦ 


T=5.0T^ ~z 






▲ 


T=10 : 


- 










































































-▲- 










■s 
















\ 
















\ 


1 1 


1 1 


1 1 1 


L^U ^ 


)-" 


-3 

10 


10' 


2 




D 








P-Pc 







Figure 3: Log-log plot of characteristic time for damage healing t, obtained by fitting equation ^ to data as shown 
in figure [21 versus p — . Data obtained using lattices of size L = 400 and different temperatures, as listed in the 
figure. The dashed line, which has slope = 1.25, is the slope obtained averaging over all measured temperatures 
and has been drawn for the sake of comparison. 
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Figure 4: Lineal-lineal plots of Dgat versus — p obtained for different values of T, as listed in the figure. 
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Figure 5: Lineal- lineal plots of the slopes obtained by fitting the data shown in figure |4] by means of equation ^ 
given by F{T), versus {T — Tc) /Tc . The straight line with slope A « 1.5 was obtained by fitting the data according 
to equation ([8|. 
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Figure 6: Scaling plot of Dgat according to equation obtained taking Fq = 0. More details in the text. 
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Figure 7: Phase diagram for the proposed model for damage spreading with a probability of damage healing. The 
full line has been drawn to guide the eyes. The dashed line corresponds to the fit of the data for T ^ 1.5Tc and 
has slope —a = — 1. 



12 



